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patients  under  a  physician's  supervision.  The  present  study  presents 
a  model  for  the  operation  of  the  Acute  Minor  Illness  Clinic  (AMIC)  part 
of  the  program  installed  at  Silas  B.  Hays  Hospital,  Fort  Ord,  California. 

The  study  draws  an  analogy  between  the  clinic  and  a  mult i- channel, 
time  dependent  queue ing  system.  An  analytic  model  is  presented  which 
uses  time  dependent  arrivals  and  server  schedules  and  determines 
expected  queue  size,  expected  delay  and  the  probability  of  delays 
exceeding  a  given  maximum  length. 
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I.  A-DSIST  PROGRAM 

A.  HISTORY,  PURPOSE,  AND  ORGANIZATION 

In  the  late  1960's  the  United  States  Army  became  aware  that  with  the 
termination  of  the  draft  a  source  of  physicians  for  the  Armed  Services 
would  be  lost.  This  loss  would  certainly  have  profound  effects  upon 
the  heal tli  services  provided  by  the  U.S.  Army,  especially  in  the  areas 
of  clinic  care.  In  1969  in  order  to  help  maintain  adequate  outpatient 
medical  care  in  the  face  of  manpower  shortages,  the  Office  of  the 
Surgeon  General  formed  the  Automated  Military  Outpatient  System  (AMOS) 
Project. 

From  the  work  of  Project  AMOS  has  come  the  AMOS  I  ST  Program.  This 
program  lias  emerged  due  to  the  ability  of  Project  AMDS  to  record  the 
logic  used  by  physicians  in  their  diagnosis  and  treatment  of  minor 
illnesses  so  that  a  nonprofessional  medical  person  (ADS  I  ST)  with  minimal 
training  could  safely  treat  certain  types  of  medical  problems.  The 
AMOSIST  would  then  become  a  provider  of  health  care,  extending  the 
capabilities  of  the  attending  physicians.  The  ADSIST  Program  is  solely 
for  ambulatory  health  care. 

The  AMOSIST  Program  is  made  up  of  two  parts  called  the  Triage  and 
the  AMIC.  In  the  Triage  part  of  the  program  an  ambulatory  patient  is 
screened  by  an  ADSIST  using  a  "Triage  Note"  (diagnostic  decision  tree 
questionnaire)  in  order  to  determine  to  which  hospital  clinic  the  patient 
should  be  sent.  The  Acute  Minor  Illness  Clinic  (AMIC)  is  staffed  by 
AMOSISTs  and  physicians  (ASSIST  MDs)  who  treat  people  screened  for 
the  AMIC  by  the  Triage.  The  "Triage  Note"  will  specif)'  whether  the 


patient  is  to  see  an  ADSIST  or  an  A*DSIST  MD,  but  a  patient  may  request 
to  sec  only  an  ASSIST  MD  if  he  desires. 

AvDSISTs  are  always  under  the  supervision  of  a  physician,  who  sees 
patients  sent  directly  from  the  Triage  Section  and  patients  referred 
from  the  various  AvDSISTs.  An  ADSIST  may  refer  to  a  physician  if  he 
does  not  feel  sure  of  his  diagnosis  or  is  faced  with  an  unfamiliar 
complaint.  The  physician  may  merely  confer  with  the  ADSIST  verbally 
or  he  may  go  to  the  ADSIST 's  room  to  see  the  patient  himself. 

In  order  to  make  the  program  successful,  the  criteria  for  AMOSISTs 
are  very  stringent.  Upon  selection,  these  personnel  are  given  two  weeks 
didactic  training  at  Fort  Sam  Houston,  Texas,  followed  by  about  ten 
weeks  on  the  job  training  as  stated  in  Reference  1.  Physicians  attend 
the  two  weeks  of  didactic  training  and  may  then  return  to  their  facilities 
to  help  ready  the  clinic. 

At  present  the  AMOSIST  Program  lias  been  adopted  in  over  20  Army 
Medical  Facilities.  Adoption  of  this  program  at  a  particular  facility 
is  decided  after  studies  have  been  done  to  ensure  the  program's  feasibility. 

B.  IMPLEMENTATION  AT  SILAS  B.  HAYS  HOSPITAL 

In  December,  1973,  the  ADSIST  Program  was  instituted  at  Silas  B. 
Hays  Hospital,  Fort  Ord,  California.  Figure  1  shows  the  physical  layout 
of  the  A IIC. 

The  Triage  is  physically  separate  from  the  AMIC  and  is  located  next 
to  the  patient  records  room.  Patients  who  are  directed  (Triaged)  to  the 
AMIC  arc  logged  in  at  the  AMIC  reception  desk  and  have  necessary  vital 
signs  taken  by  a  Nursing  Assistant.  The  patient's  folder  is  then  placed 
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in  either  a  physician  or  ATDSIST  rack  depending  on  the  "Triage  Mote" 
attached.  'Hie  patient  is  seen  when  a  server  becomes  available.  Order 
of  service  is  first  in,  first  out. 

During  a  treatment,  a  chaperone  may  be  necessary.  Nursing  Assistants 
or  even  ADSISTs  are  used  for  this.  If  referral  to  a  physician  is  nec- 
essary, the  ATDSIST  will  look  for  an  available  physician.  He  may  have 
to  wait  if  all  physicians  are  busy.  When  a  physician  is  located,  either 
a  verbal  referral  will  occur  or  the  physician  will  go  with  the  AMOS 1ST 
to  examine  the  patient.  In  either  case  the  effect  is  the  same,  namely 
that  two  persons  are  occupied  with  one  patient. 

A  small  percentage  of  patients  come  to  the  AMI C  just  for  prescription 
refills.  Their  records  are  placed  in  a  separate  rack.  Refill  patients 
are  given  a  priority  by  the  physicians  and  wait  in  a  separate  area  from 
the  other  patients.  Their  service  times  arc  normally  very  much  shorter 
than  those  for  other  patients. 

Another  small  percentage  of  patients  are  sent  from  the  AMIC  to  the 
X-RAY  Department  or  to  the  Laboratory.  Upon  their  return  to  the  AMIC, 
they  are  given  priority  by  the  AMOSIST  or  physician  who  originally 
saw  them. 

Figure  2  presents  a  block  diagram  of  the  system  just  described.  It 
can  be  determined  at  this  point  that,  with  the  exception  of  referrals 
from  the  AYQSISTs  to  physicians,  the  system  can  be  viewed  as  two  separate 
multi-channel  queues.  The  analytical  model  will  be  discussed  in  another 
Section. 
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II.  AMIC  DATA 

In  order  to  develop  an  analytical  model  capable  of  analyzing  and 
predicting  the  system  behavior  of  the  AMIC  at  Silas  B.Hays  Hospital  it  was 
necessary  to  establish  a  data  base.  To  keep  the  modeling  of  the  AMIC 
in  a  workable  form,  all  statistics  related  to  prescription  refills  and 
patients  sent  from  the  AMIC  to  X-RAY  or  IAB  were  disregarded.  Also, 
complete  statistics  were  confined  to  weekdays  when  patient  load  was  the 
heaviest. 

A.  PERSONNEL  STAFFING  AND  SCHEDULING 

During  this  study,  there  were  four  physicians  available,  one  being 
responsible  for  the  clinic's  medical  administration  as  well  as  seeing 
patients.  There  were  11  AMOSISTs  available  who  either  saw  patients  or 
assisted  in  the  Triage.  An  additional  non-medical  enlisted  person  was 
the  NCOIC.  The  support  staff  was  made  up  of  four  Nursing  Assistants 
and  three  receptionists. 

The  AMIC  was  open  to  receive  patients  from  0745  to  2330  hours, 
with  active  duty  personnel  having  priority  prior  to  1000  hours.  The  day 
shift  of  personnel  work  until  1630  while  the  night  shift  worked  from 
1530  to  2330  hours.  During  weekdays,  three  AMOS 1ST  MDs  worked  on  the 
day  shift  and  one  at  night.  Since  the  AMIC  was  tasked  with  sending  a 
physician  to  the  stockade  each  morning,  this  meant  that  at  most  two 
physicians  would  be  available  during  the  mornings. 

On  Monday,  the  heaviest  day,  seven  to  eight  AMOSISTs  would  be 
scheduled  to  work  the  day  shift  and  on  the  other  weekdays,  five  to  six. 
At  night,  two  AMOSISTs  were  scheduled  on  all  weekdays. 


13 


B.  PATIENT  ARRIVAL  PATTERN 

In  order  to  get  an  adequate  estimate  of  the  patient  arrival  pattern 
in  the  AMIC,  log-in-sheets  for  February  to  May  1974  were  reviewed.  Each 
month's  arrivals  were  recorded  for  each  hour  of  the  day  and  for  each 
day  of  the  week.  These  arrival  figures  were  averaged  together  to  give 
a  standard  week.  Holidays  and  days  affected  by  holidays  were  left  out 
of  the  calculations  since  they  were  unusual  in  value  and  not  representative 
of  an  average  day.  For  periods  that  did  not  include  a  full  hour  period, 
the  figures  were  adjusted  to  an  hourly  figure  in  order  to  make  them 
comparative  for  statistical  purposes .  Table  I  summarizes  the  arrival 
figures  for  an  average  week.  Each  day's  average  was  based  on  17  weeks 
of  data. 

Each  day  of  the  week  exhibits  approximately  the  same  shape  curve  for 
arrival  patterns  over  the  length  of  the  day  and  differs  only  in  the 
average  daily  total  number  of  arrivals.  This  allowed  arrival  figures  to 
be  pooled  together  giving  one  average  weekday  arrival  pattern  as  depicted 
in  Figure  3.  Different  patient  volumes  could  then  be  achieved  by  merely 
moving  the  entire  curve  up  or  down  to  produce  any  desired  patient  volume. 

C.  AMOSIST  MD  SERVICE  TIMES 

There  was  no  procedure  within  the  AMIC  to  recover  the  length  of  time 
that  an  A'DSIST  MD  spends  with  a  patient  or  with  an  AMOSIST  on  a  referral. 
To  estimate  these  times,  approximately  75  actual  observations  of  patient 
encounters  with  A'OSIST  MDs  were  made  during  randomly  selected  days  over 
a  three  week  period.  The  results  of   these  observations  showed  that  an 
AMOSIST  MD  spends  approximately  13.55  minutes  with  his  own  patients  and 
about  4.34  minutes  of  his  time  with  an  AMOSIST  on  a  referral. 
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0745  -  0800* 

0800  -  0900 

0900  -  1000 

1000  -  1100 

1100  -  1200 

1200  -  1500 

1300  -  1400 

1400  -  1500 

1500  -  1600 

1600  -  1700 

1700  -  1800 

1800  -  1900 

1900  -  2000 

2000  -  2100 

2100  -  2200 

2200  -  2500 

2300  -  2530* 


TABLE  I 
ARRIVALS  PER  HOUR  FOR  AVERAGE  WEEK 

MON.  T!JES.  V.ED.  TILT..  FRI . 


18.48 

15.24 

10.36 

9.88 

12.83 

24 .  69 

23.50 

25.12 

22.00 

20.67 

21.75 

20.81 

17.94 

16.65 

14.94 

17.50 

13.56 

16.47 

13.76 

12.56 

10.25 

9.00 

8.41 

S.71 

9.67 

7.88 

7.56 

6.12 

5.94 

6.72 

11.31 

13.50 

12.47 

12.65 

11.06 

10.19 

11.13 

9.71 

11.82 

9.44 

9.06 

8.31 

7.65 

7.82 

8.50 

8.06 

5.69 

6.24 

6.82 

6.11 

7.31 

6.56 

5.12 

4.76 

4.33 

5.44 

5.94 

6.71 

6.47 

4.39 

6.25 

6.69 

5.71 

5.65 

4.67 

5.13 

C      70 
O  .  .Jo 

4.47 

4.76 

4.33 

3.56 

3.88 

3.06 

3.24 

2.83 

2.00 

2.50 

2.55 

3.11 

3.22 

.50 

.76 

.88 

.88 

1.10 

*     ADJUSTED  TO  HOURLY  AVERAGE 
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D.  ASSIST  SERVICE  TIMES 

The  service  times  for  each  Af-DSIST  were  broken  down  into  referral 
and  non-referral  categories  just  as  for  the  physicians.  The  actual 
times  spent  with  patients  for  each  individual  A'OSIST  were  taken  from 
the  Data  Collection  Sheets  (computer  generated  questionnaire  used  by 
A-DSIST  to  diagnose  a  patient's  complaint)  which  are  kept  on  file  for 
each  A^DSIST.  These  sheets  also  show  by  a  physician's  signature  or  an 
appropriate  box  being  checked  when  a  referral  was  made.  After  analysis 
of  these  actual  patient  encounter  times,  an  average  service  time  (time 
from  start  of  one  patient  to  start  of  next  patient)  was  derived  for 
each  ADSIST.  In  all  but  one  case,  at  least  40  encounters  were  used  to 
compute  average  "service  times"  for  an  ADSIST. 

The  results  of  this  data  indicated  that  the  average  service  time  for 
a  non-referred  patient  was  17.96  minutes  while  the  referred  patient's 
service  time  was  23.85  minutes.  The  actual  patient  encounter  times  for 
the  ADSISTs  compared  quite  favorably  with  results  published  in  a  recent 
ADS  Project  Progress  Report  (Reference  2).  In  addition,  an  average 
referral  rate  for  an  AFDSIST  was  found  to  be  30.93°6.  The  above  estimates 
were  based  on  about  400  AFDSIST' s  encounters. 

E.  '  01I1ER  STATISTICS 

Since  the  AMIC  is  a  system  that  generates  patients  to  two  different 
types  of  servers,  it  was  necessary  to  find  the  percentage  of  people  who 
were  either  triaged  to  sec  an  AFDSIST  MD  or  requested  to  see  an  ASSIST 
MD.  This  type  of  data  was  not  readily  available  from  the  records  of 
the  AMIC.  For  selected  weeks  the  Data  Collection  Sheets  for  all  AMOSISTs 
were  totaled  for  certain  days.  This  figure  showed  the  total  patients 
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triagcd  to  an  ADSIST  who  actually  saw  an  AMOSIST  on  a  given  day.  Tliis 
was  subtracted  from  the  total  patients  entering  the  AMIC,  as  shown  on 
the  log-in-sheets,  to  give  the  number  of  patients  seeing  an  A*DSIST  MD 
directly.  This  figure  showed  that  approximately  55%  of  the  patients 
were  sent  directly  to  the  ADSIST  MD  versus  45?,  for  the  ANDSIST.  This 
estimate  originated  from  approximately  1000  cases. 

Since  it  was  decided  to  start  the  analytical  model  at  0800  in  the 
morning,  an  estimate  of  the  average  number  of  people  in  the  waiting 
area  prior  to  the  first  patient  being  seen  was  required.  The  arrival 
information  previously  gathered  for  February  to  May,  1974,  was  used  to 
find  an  estimate  of  11  people  waiting  to  see  a  server  prior  to  the  first 
sendee  time. 

F.  INTERPRETATION  OF  DATA 

It  would  be  inaccurate  to  say  at  this  point  that  the  AMIC  at  Silas  B. 
Hays  Hospital  can  be  fully  represented  with  the  aforementioned  set  of 
numbers.  Any  system  that  is  composed  of  people  and  is  dependent  upon 
their  interactions  can  probably  never  be  completely  and  accurately 
modeled. 

The  data  gathered  on  the  system  were  required  by  the  analytical 
approach  to  the  AMIC  model.  Also,  the  analytical  model  is  a  long  run 
representation  and  average  figures  should  not  be  expected  to  predict 
exactly  what  the  situation  will  be  like  tomorrow.  Some  variation  will 
almost  always  occur. 

Each  system  has  its  own  inherent  characteristics  that  may  not  be 
readily  discernible  at  first  glance.  The  AMIC  at  Silas  B.  Hays  is  no 
exception.  The  stockade  commitment  in  the  mornings  requiring  a  physician 
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from  the  AMIC,  the  priority  until  1000  for  active  duty  personnel  giving 
an  uneven  arrival  pattern,  and  the  makeup  of  the  population  that  the  AMIC 
serves,  arc  things  that  may  not  be  common  to  any  other  AMIC.  Each 
system  must  be  judged  on  its  own  performance  in  the  face  of  its  individual 
problems  and  not  be  compared  using  "raw  data"  to  other  systems  which 
outwardly  appear  the  same. 

The  analytic  model  output  to  be  presented  pertain  to  Silas  B.  Hays 
Hospital  as  it  was  operating  during  the  time  of  this  study  and  uses 
data  taken  at  a  particular  period  of  time.  No  doubt  as  the  system 
changes  in  time  certain  statistics  may  need  revision.  The  greatest  use 
of  this  study  will  be  the  chance  to  see  the  results  of  changes,  either 
proposed  or  actual,  and  predict  their  effect  upon  the  entire  system. 

Table  II  summarizes  the  data  that  has  been  presented  in  this 
Section  and  which  will  be  used  as  the  data  base  for  the  analytic  model. 
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TABLE  II 


DATA  SUMMARY 


AMOSIST  Ml): 

MEAN  SERVICE   (DIRECT) 
MEAN  SERVICE   (REFERRAL) 

AIvDSIST: 

MEAN  SERVICE  (NON -REFERRAL) 
MEAN  SERVICE  (REFERRAL) 
REFERRAL  RATI: 

PROPORTION  OF  DIRECT  :iD  PATIENTS 
MEAN  ARRIVALS  PRIOR  TO  FIRST  SERVICE 
AVERAGE  ARRIVALS  PER  WEEKDAY: 

MON.         TUES.         WED. 

156         149         145 


13.55  min. 
4.34  min. 


17.96  min. 

23.83  min. 

30.93% 

55.00% 

11 

TIIURS. 

FRI 

140 

127 

STAFFING: 

A^DSIST  MD 
ASSIST 
NURSING  ASST. 
NCOIC 
RECEPTIONIST 


4 
11 

4 
1 
3 
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III.  Till!  AMIC  AS  A  QUHJEING  SYSTB1 

During  the  course  of  data  gathering  and  research,  the  similarities 
of  the  Acute  Minor  Illness  Clinic  to  a  multi- channel  queue ing  system 
became  apparent.  In  the  actual  system  the  arrivals  follow  an  inhomo- 
geneous  pattern  throughout  the  course  of  any  day.  Once  in  the  system 
a  patient  is  directed  to  an  AMOSIST  or  an  ADSIST  MD  and  then  waits 
for  the  first  available  server.  l\hen  a  server  is  free,  the  earliest 
arrival  waiting  is  treated  and  released.  The  major  complication  in 
modeling  the  actual  system  is  the  interference  in  the  AMOSIST  MD  queue 
caused  by  referrals  from  AMOSISTs.  However,  if  the  two  queues  (AMOSIST 
and  AMOSIST  MD)  can  be  isolated,  two  independent  queueing  problems  can 
be  formulated. 

The  Af  JOS  I  ST  referrals  cause  an  interference  in  the  AMOSIST  MD  queue 
because  the  referred  patient  occupies  an  AMOSIST  MD  with  an  additional 
service.  This  however,  is  the  only  effect  because  an  AMOSIST  remains 
occupied  while  the  referral  takes  place.  Thus  the  problem  can  be 
simplified  by  determining  an  appropriate  augmentation  to  the  AMOSIST  MD 
arrival  stream  to  handle  those  patients  referred  by  AMOSISTs.  These 
additional  arrivals  can  be  considered  as  "virtual"  arrivals  -  i.e. 
persons  who  will  eventually  occupy  the  AMOSIST  MD's  time  but  who  initially 
arrive  to  see  the  AMOSIST. 

The  referral  rate  from  A-DSIST  to  ADSIST  MD  is  taken  as  a  percentage 
of  the  ADSIST  arrival  stream.  The  number  of  arrivals  thus  computed  was 
then  added  to  the  ADSIST  MD  arrival  stream.  The  original  clinic  system 
can  then  be  viewed  as  two  distinct  multi-channel  queueing  systems. 


In  order  to  use  the  data  presented  in  Section  II  as  an  input  to  the 
analytic  model,  it  was  necessary  to  calculate  certain  new  statistics. 
The  hourly  arrivals  as  well  as  a  single  average  service  time  for  an 
AM0S1ST  and  an  A'DSIST  MD  were  computed. 

For  a  given  hourly  arrival  figure,  the  product  of  this  number  and 
the  average  percentage  seeing  an  A^DSIST  gives  the  hourly  arrivals  for 
the  AMOSIST.  The  hourly  arrivals  for  an  A-DSIST  MD  is  the  product  of 
total  hourly  arrivals  and  the  percentage  seeing  the  physician  directly 
added,  to  the  product  of  the  referral  rate  and  the  hourly  arrivals  for 
an  ADS  I  ST. 

The  analytic  model  recognizes  only  one  average  service  time  for  an 
ADSIST  or  an  AMOSIST  MD.  For  both  servers,  Section  II  pives  two 
different  types  of  service  times.  For  the  A'DSIST,  it  was  referral 
and  non-referral  service  times.  'lrhc   AIOSIST  MD  had  either  direct  patient 
or  referral  service  times.  The  percentage  of  each  type  of  service  given 
by  an  AMOSIST  or  AIOSIST  MD  was  computed.  The  product  of  this  percentage 
and  the  appropriate  average  service  time  added  over  types  of  service 
gives  an  effective  average  service  time  for  the  A'DSIST  and  ADSIST  MD. 
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IV.  TIDi  MODEL 

A.  GENERAL  DESCRIPTION 

In  order  to  analyze  the  two  queucing  systems,  a  mathematical  model 
was  developed.  Using  the  inhomogeneous  arrival  stream  and  server  schedules 
as  input,  the  model  computes  the  probability  distribution  of  system  size 
throughout  the  day.  The  model  then  computes  various  measures  of  system 
performance  at  specified  times  of  the  day.  The  model  treats  the  two 
queues  as  separate  systems  and  during  each  day  solves  first  the  A*DSIST 
system  and  then  the  A'DSIST  MD  system. 

The  significant  problem  in  the  formulation  of  the  model  was  that 
all  rates  and  schedules  of  servers  were  functions  of  time.  There  is  no 
known  analytical  method  of  solution  for  the  time  dependent  multi- server 
queue ing  system.  Katadare  and  Kaufman  (Reference  3)  have  advanced  a 
Method  of  Successive  Approximate  Transient  Solutions  (MDSATS)  for 
numerical  analysis  of  a  single  channel  time  dependent  system  and  much  of 
their  method  has  logical  extensions  to  the  multi-channel  system.  The 
approach  calls  for  discretizing  the  time  axis  into  segments  over  which 
the  parameters  remain  constant.  Then  a  transient  solution  is  computed. 
for  each  segment  successively.  The  method  was  expanded  to  multiple 
channel  queues  and  made  efficient  for  real  time  computations  for  its 
application  to  this  problem.  The  details  are  explained  more  fully  below. 

B.  MODEL  ASSUMPTIONS 

The  model  makes  several  key  simplifying  assumptions.  These  will  be 
listed  and  explained  briefly  here  and  referred  to  in  subsequent  sections. 
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(Al)  Arrivals  form  an  inhomogencous  Poisson  Process. 

(A2)  Services  are  exponentially  distributed  and  independent  of 
the  arrival  stream. 

(A3)  At  most  three  arrivals  and  three  departures  can  occur  in  any 
micro  time  step.  A  micro  time  step  is  a  sub-interval  of  the  periods 
in  which  parameters  remain  constant;  it  is  used  in  the  computation  of 
state  probabilities. 

(A4)  No  patient  can  both  arrive  and  depart  within  the  same  micro 
time  step.  Patients  leaving  during  a  micro  time  step  are  from  those 
present  at  the  beginning  of  the  step.  This  assumption  has  the  effect 
•  of  simplifying  the  calculations  somewhat. 

C.  r>DDPL  OPERATION 

In  order  to  adequately  describe  the  operation  of  the  model,  a  general 
description  of  the  computations  for  solving  a  day's  clinic  operation  is 
presented  here.  Mathematical  amplifications  are  presented  with  the 
appropriate  steps  where  necessary.  Output  and  interpretation  are  discussed 
later. 

1.   Initial  State  Probability  Vector 

The  probability  distribution  of  system  size  at  the  day's  beginning 

is  Poisson  with  a  mean  given  by  the  user.  Specifically, 

P  (0)  =  P  ,c     n  =  0  1  ?     \T 
nl  J       ~rrt >  ,l   u>J-»->  •  •  •  >iN 

where  Pn(t)  is  defined  to  be  the  probability  that  system  size  is  n  at 
time  t  and  p  is  the  average  number  of  early  arrivals  for  the  given  day. 
To  ensure  that  no  probability  mass  was  lost  by  limiting  system  size  to 
N,  the  model  places  all  remaining  mass  in  the  probability  of  largest 
system  size.  That  is, 
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N-l 

PN(0)  =i.o-2  Pn(0). 
H  n=0 

In  all  applications  this  value  was  zero.  System  size  (N)  used  in  the 

model  was  15  for  computational  efficiency.  This  restriction  is  not 
severe  because  in  light  of  the  two  queue  concept,  it  represents  an  over- 
all system  size  of  30. 

2.  Time  Step  Computations 

The  model  next  determines  the  first  period  of  time  during  which 

all  parameters  (arrival  rate  and  number  of  servers)  arc  constant.   It 

should  be  noted  that  both  arrival  rates  and  server  schedules  are  specified 

as  step  functions.  Within  this  interval,  subsequently  referred  to  as  the 

macro  interval,  time  dependence  is  removed.  Using  (A3)  and  (A4)  this 

macro  interval  is  subdivided  into  micro  intervals  such  that  the  assumptions 

are  valid.  Specifically, 

f  _  M._  r1.25    1.25  . 

where  &   is  the  length  of  the  micro  interval,  A  is  the  appropriate 
arrival  rate,  c  is  the  number  of  servers  and/*  is  the  service  rate  per 
server.  The  length  of  the  macro  interval  is  then  divided  by  o  to 
determine  the  number  of  "steps"  which  will  be  necessary  to  compute  the 
state  probability  vector  at  the  end  of  the  period.  The  constant  (1.25) 
used  in  the  6    computation  was  experimentally  determined  as  that  figure 
which  gave  the  best  approximation  to  a  correct  steady  state  solution  for 
a  time  independent  system. 

3.  Transition  Probability  Vatrix 

In  order  to  determine  the  probability  of  entering  a  given  state 

from  an  initial  state  in  a  o  step,  a  matrix  of  transition  probabilities 

[Q]  is  next  determined. .  The  entries  in  this  matrix  arc  defined  as 

Ok  n^)  =  ^r  (system  size  at  t  is  n  given  system  size  at 
'         t-S     was  k) 
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Further  define 

Q-,  •  •  (t)  =  Pr  (i  arrivals  and  j  departures  in  6    given 
'  '•'        system  size  at  t-«$   was  k). 

Then  we  have  approximately  that 

Qk>nw  *  TL     %i>jm 

l  *  j 

j  *<  3 
k+i-j=n 

To  determine  the  Q,   two  more  definitions  will  be  necessary: 

K  ,11 

X(i)  =  Pr  (i  arrivals  in  a  $    step)  i  =  0,1,2,3 

Y(j  ,1c)  =  Pr  (j  departures  in  a  S     step  given  k  present 

at  the  beginning  of  the  6    step)  j  =  0,1,2,3. 

It  is  clear  that  the  probability  of  an  arrival  is  independent  of  system 

state  and/or  number  of  servers.  In  fact,  by  (Al) ,  the  probability  of 

i  arrivals  in  6    is 

X(i)  =  (M)  e    /V.      ,  i  =  0,1,2,3 
By  (A3)  the  number  of  departures  in  6   is  limited  to  three.  This 
restriction  defines  three  distinct  cases  for  the  probability  of  a 
departure  based  on  system  size  and  number  of  servers.  These  cases  are 
first  discussed  and  then  the  combination  of  arrival  probabilities  and 
departure  probabilities  to  determine  the  Qi,  n(t)  entries  is  discussed. 
a.  System  Size  Less  Than  c 

An  entering  system  size  less  than  the  number  of  servers 
implies  that  not  all  servers  arc  occupied  and  a  binomial  probability  of 
departure  exists.  Specifically, 

YU.10  ={■)  Cl.O-e*"  )hcM)k~i    ,  3  -   0,1,2,3  and  k>,  j . 
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b.  System  Size  Greater  Than  c+2 

When  the  beginning  system  size  exceeds  the  number  of  servers 
by  two  or  more,  the  system  is  saturated  and  no  server  becomes  idle 
during  the  period  for  cases  of  at  most  three  departures.  The  probability 
of  a  departure  is  therefore  k- independent.  Specifically, 

i  -Q.JUO 

Y(j,k)  =  (c/u6  )Je     /j!   ,  j  =  0,1,2,3  and  k>c+2. 

c.  System  Size  Between  c  and  c+2 

When  the  beginning  system  size  is  such  that  all  servers  are 
occupied  but  also  might  become  idle  at  some  point  during  the  time 
interval,  the  exponential  assumption  (A2)  must  be  interpreted  properly 
to  account  for  the  possibility  of  idle  servers.   If,  in  fact,  the  servers 
become  idle  during  &       ,  the  distribution  of  inter- departure  time  is 
2-Erlang  followed  by  3-Erlang,  etc.  depending  on  the  specific  number  of 
servers  becoming  idle.  Looking  more  closely  at  the  specific  cases 
involved,  it  can  be  shown  that: 

for  k  =  c+1,  Y(0)  =  c"C/l 

Y(l)  =  (c/.kS)  c'CM* 

Y(2)  =  c  (e  v   ^  )(l-e    -/i<Jc   ) 

Y(3)  =  1.0  -  £  Y(j) 
j=0 


and  for  k  =  c+2, 

C/A</ 


Y(0)  =  e 

Y(l)  =  (c/t<&)   c"C/Ucf 

Y(2)  =  {Qjml)2c'ZM6/2\ 

Y(3)  -  1.0  -  «£   Y(j) 
P 
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Finally,  the  entries  in  the  transition  probability  matrix  may  be 
determined.  By  (A4)  arrivals  and  departures  are  independent  events  and 
thus  Q,  .  .  =  X(i)  Y(j,k).   It  should  be  emphasized  that  the  effects  of 
k  (entering  state)  are  in  the  Y(j)  term. 

A  necessary  property  of  transition  probability  mn.tr ices  in  general 
is  that  the  row  elements  sum  to  1.0.  That  is,  given  any  beginning  state, 
the  system  must  be  in  some  state  at  the  end  of  the  interval.  By  (A3) 
the  model  discards  any  tail  probability  beyond  three  arrivals  and/or 
three  departures.  To  correct  for  this  assumption  and  to  ensure  that 
correct  row  sums  are  obtained,  the  convention  of  placing  all  remaining 
probability  mass  in  the  X(3)  and  Y(3,k)  terms  was  adopted,  except  in 
those  cases  where  the  entering  state  was  less  than  three.  In  those  cases 
the  remaining  probability  was  placed  in  the  highest  possible  number  of 
departures  consistent  with  entering  system  size.   (Recall  (A4)). 

4.  The  Interval  Computation 

Having  determined  the  6  step  transition  probability  matrix,  it 
is  necessary  to  show  the  computation  used  to  determine  follow-on  state 
probabilities.   It  can  be  readily  seen  that: 

for  P(t)  =  (P  ft),  P  (t),...,P  (t)), 
P(t+*  )  -  P(t)[Q] 
P(t+2  4  )  =  P(t+<£  )[Q]  =  P(t)[Q][Q]  =  P(t)[Q]2 


P(H«J  )  =  PWfQ]1  , 
where  i  is  the  previously  computed  number  of  6    steps  Contained  in  the 
macro  interval  under  consideration.   (See  IV. C. 2  above).  Within  the 
model,  i  is  taken  to  be  a  power  of  two  for  'computational  efficiency. 
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5.   Iterations 

The  7'csults  of  the  preceding  computations  are  state  probability 
vectors  for  the  end  of  each  macro  interval.  This  state  probability 
vector  replaces  the  initial  probability  vector  (1)  and  the  confutations 
(2  through  4)  are  repeated  iterativcly  until  a  predetermined  end  of  day 
is  reached.  The  model  then  returns  to  the  second  arrival  rate  and  server 
schedule  (/CDSIST  MD)  and  repeats  the  entire  procedure  until  the  pre- 
determined end  of  day  is  again  readied.  This  then  represents  the  em\   of 
one  day's  clinic  operations. 

D.  STATISTICS 

1.   Identified 

It  was  initially  recognized  that  probabilities  of  system  size 
and  delay  would  be  necessary  if  the  model  were  to  be  of  any  practical 
value  to  the  user.  After  considering  the  myriad  of  available  interpretive 
tools,  the  list  was  ultimately  narrowed  to  eight  distinct  measures  of 
the  system  which  are  determined  at  each  time  the  user  desires  a  report. 
These  statistics  are  defined  and  discussed  separately  in  the  following 
sections. 

a.  Expected  System  Size  is  computed  as 

£  nPn(t). 

n=0 

b.  Standard  Deviation  of  System  Size  is  computed  as 


n"pn(t0  -  (2;nPn(t))   . 
itfD  n=0 

c.  The  probability  a  server  is  busy  at  t  is 

Q  ■  (Cr  »l>  (O+c  r*   P  (t))/c. 
n=0         n=c   n 
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It  should  be  noted  that  this  probability  is  analagous  to  the  standard 
utilization  determination  of  qucucing  theory  and  is  hereafter  referred 
to  as  utilization  (or*,)« 

d.  The  expected  size  of  the  queue  at  t  is 

A  nP-(t)  "  cd  •*»«  " c  t  p»w  •     . 

e.  The  standard  deviation  of  queue  size  at  t  is 


v< 


JL       2 

JZ      (n-c)^Pn(t). 

n=c 

f .  The  probability  an  arrival  at  t  encounters  no  delay  is 

Si   P(t)   . 
n=TT   n^  J 

g.  The  probability  an  arrival  at  t  encounters  delay  in  excess 
of  d  minutes  is 


p„(t). 


n=c  I  k=0      *.         ) 

h.  The  expected  delay  for  an  arrival  at  t  is 

N          c-1          N  N 

71     nPn(t)  -  Yl     "PnW  "  cIPnW  +  H  Pn(t) 

n=D         n=0          n=c  n=c 
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2.     Interpretation 


It  is  worthwhile  to  note  that  because  of  the  augmentation  of  the 
A^OSIST  J  ID  arrival  rate  it  is  necessary  to  adjust  certain  of  the  above 
statistics  to  remove  the  effect  of  the  "virtual"  arrivals.  Statistics 
a,b,d  and  c  are  therefore  interpreted  as  follows:  a  and  d:  the 
expectation  of  system  and  queue  size  are  adjusted  by  a  factor  of  p  where 
p  represents  the  percentage  of  the  arrival  stream  which  are  actual 
A>DSIST  MD  arrivals;  b  and  e:  the  standard  deviations  are  computed  as 


\£ 


;(X)  +  p2V(X) 
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where  p  is  as  above,  q  is  (1-p)  and  X  represents  the  computed  and 
unadjusted  figures  for  b  and  e. 

E.  IBM  360/67  APPLICATION 

1.  General 

The  above  model  has  been  programmed  in  G-Level  Standard  Fortran 
for  implementation  on  the  IBM  360/67  under  OS  and  for  the  CP/CMS  time 
sharing  system  at  Naval  Postgraduate  School.  The  general  flow  chart  of 
the  model  is  included  as  Figure  4.  The  model  successfully  reflects  the 
actual  system  operation  of  the  AMIC  at  Silas  B.  Hays  Hospital.  It  has 
been  incorporated  in  an  interactive  time  sharing  mode  and  through  the 
use  of  portable  terminal  equipment  has  been  demonstrated  to  hospital  and 
clinic  administrators. 

2.  Input 

The  model  is  sufficiently  general  to  allow  user  control  of  the 
following  input  parameters: 

a.  Volume. 

The  user  may  specify  any  number  of  arrivals  desired  in  the 
clinic  for  a  given  day.  The  model  interprets  this  volume  as  a  scale 
factor  for  the  standard  arrival  rate  step  function. 

b.  Service  Time 

The  user  may  vary  the  mean  service  time  for  cither  type  of 
server  and  also  for  referral  and  non-referral  service. 

c.  Schedules 

The  user  may  specify  schedules  for  each  type  of  server  and 
may  have  the  number  of  each  type  server  on  duty  change  as  many  as  16 
times  throughout  the  day. 
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d.  Percent  of  Arrival  Stream 

The  user  may  specify  that  percentage  of  the  arrival  stream 
which  sees  the  doctor  directly. 

e.  Referral  Rate 

The  user  may  specify  that  percentage  of  patients  which  will 
ultimately  see  a  doctor  through  referral  from  an  ASSIST. 

f.  Early  Arrivals 

The  user  may  specify  the  number  of  arrivals  before  the  system 
opens . 

g.  Reports  Desired 

The  user  may  specify  up  to  16  times  throughout  the  day  at 
which  he  would  like  statistics  computed  and  displayed. 

h.  Delay  Cutoff 

The  user  may  specify  the  number  of  minutes  for  which  he 
would  like  the  probability  of  delay  in  excess  of  that  figure  computed, 
(the  "d"  in  IV.D.l.g  above). 
3.  Output 

Although  eight  statistics  are  generated  within  the  model,  it  is 
apparent  that  not  all  of  these  have  meaning  for  the  non  statistically 
oriented  user.  For  this  reason  only  five  of  the  eight  are  normally 
displayed  in  interactive  user  sessions.  These  five  are  time  of  the 
report,  utilization,  expected  size  of  the  queue,  expected  delay  encountered 
and  the  probability  of  delay  in  excess  of  the  user  supplied  delay  cutoff. 
11  ic  entire  list  of  statistics  is  printed  off-line. 

Tables  3  and  4  show  the  statistics  generated  for  a  typical 
weekday  for  A»-DSIST  and  AJDSIST  MD  respectively.  The  server  schedules 
are  also  shown.  The  normal  method  of  gathering  statistics  is  to  compute 
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at  one  minute  past  cacli  hour  to  capture  conditions  which  will  be  in 
force  for  the  coming  period.  That  is,  because  most  changes  in  schedule 
and  arrival  rate  step  functions  occur  on  the  hour,  the  statistics  are 
taken  during  the  following  period  so  that  the  conditions  which  will 
attempt  to  clear  the  system  will  be  in  force. 
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V.  CONCLUSIONS 

A.  VALIDATION  AM)  SENSITIVITY 

The  computer  model  was  validated  by  forcing  "steady  state"  solutions 
to  be  achieved.  This  can  be  easily  done  by  assuming  a  constant  mean 
arrival  rate,  a  constant  number  of  servers  and  a  constant  mean  service 
rate.   In  so  doing,  it  became  apparent  that  the  model  is  sensitive  to 
several  parameters. 

1.  Early  Arrivals 

The  effect  of  arrivals  before  the  system  is  open  for  business 
cannot  be  overemphasized.  The  effect  of  the  early  arrivals  is  intuitive 
in  that  if  the  system  begins  congested,  then  it  must  spend  a  good  deal 
of  time  alleviating  this  congestion.  Arrivals  during  the  initial  periods 
will  normally  experience  much  longer  delays  than  if  no  patients  arrived 
early.  Within  the  model,  early  arrivals  is  an  input  parameter. 

2.  A>PSIST  MP  Schedule 

As  the  output  tables  indicate,  during  normal  operation  the  major 
system  congestion  occurs  in  the  A^DSIST  MD  queue.  This  is  attributable 
to  the  percentage  of  the  arrival  stream  which  is  directed  to  the  AMOSIST  MD 
and  to  the  augmentation  of  this  arrival  stream  by  referrals  from  the 
ASSIST  as  well  as  the  number  of  AMOSIST  MDs  available.  This  model  is 
very  sensitive  to  the  7VDSIST  MD  schedule  and  a  large  expected  delay  can 
be  seen  when  few  servers  are  scheduled  during  peak  periods.  In  actual 
clinic  operation  the  AMOSIST  MD  arrival  rate  is  affected  by  those  patients 
who  request  to  see  an  A'DSIST  'ID  and  also  by  those  A*DSISTs  who  refer  a 
higher  percentage  of  their  patients.  The  percentage  of  arrivals  going 
directly  to  the  ADSIST  MD  and  the  referral  rate  from  AMOSISTs  to  A'DSIST 
MDs  are  inputs  to  the  model. 
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B.  SIP-MOT 

v  The  Acute  Minor  Illness  Clinic  at  Silas  B.  Hays  Hospital  can  be 
successfully  modelled  as  a  multi- channel,  time  dependent  queueing  system 
and  approximate  solutions  to  desired  queue  and  delay  statistics  can  be 
found  using  numerical  methods. 

The  computerized  version  of  the  model  is  capable  of  handling  much 
larger  systems  and  could  be  expanded  to  handle  more  types  of  servers. 
As  previously  mentioned,  the  present  version  with  system  size  limited 
to  15  is  used  interactively  with  the  user  supplying  input  and  seeing 
output  statistics  within  30  seconds. 

A  general  method  of  solution  for  multi -channel,  time  dependent 
queueing  systems  has  been  presented.  This  model  could  find  ready  appli- 
cation to  any  number  of  similar  queueing  systems  presently  unsolved. 
The  procedure  presented  is  sufficiently  general  so  that  a  wide  variety 
of,  and  variance  in,  input  parameters  may  be  easily  handled  thereby 
enhancing  the  model's  applicability. 
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